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Abstract 

I prove an identity between the first kind and the third kind complete 
elhptic integrals with the following form: 
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This relation can be applied to eliminate the complete elliptic integral of 
the third kind from the analytic solutions of the imaginary part of two-loop 
sunset diagrams in the equal mass case. 

The validity of this relation in the complex domain is also briefly dis- 
cussed. 
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Elliptic integrals, as one of the most important classes of nonelementary 
functions, have found numerous applications in many branches of engineer- 
ing and physics. Having been comprehensively studied by many eminent 
mathematicians over centuries, a vast amount of knowledge on the elliptic 
integrals, as well as their close cousin, elliptic functions, has been accumu- 
lated IHElElll]. 

The simplest and particularly important class of elliptic integrals is the 
complete elliptic integrals, totally of three kinds. The complete elliptic inte- 
gral of the third kind, 11, being the most complicated one, can be expressed 
in terms of the complete elliptic integral of the first kind, K, plus elementary 
functions and Heuman's Lambda function or Jacobi's Zeta function. The last 
two transcendental functions can in turn be expressible from the incomplete 
elliptic integrals of the first and second kinds. 

The aim of this note is to establish a simple relation between the complete 
elliptic integrals of the first and the third kind, for some specific texture 
of arguments of course. That is, in this special case, the complete elliptic 
integral of the third kind can be transformed to the complete elliptic integral 
of the first kind, plus elementary function at most, without resort to any 
other nonelementary function. 

There are a variety of conventions adopted in literature in defining the 
elliptic integrals. In this work I find it convenient to use the same definitions 
as taken by [21 |3] for the first, second and third kind of complete elliptic 
integrals: 

K(m) = [ dt , ^ = - F (-,-;l;m] , 
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where F and Fi denote Gaussian hypergeometric function, and Appell hy- 
pergeometric function of two variables, respectively. In most practical ap- 
plications, the parameter m and the characteristic n are restricted to be 
less than 1. However, it is worth emphasizing when these arguments exceed 
1 or even are off the real axis, these elliptic integrals are still well defined 
mathematically, though become complex valued in general. 
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The main result states as follows. For Vx G i?, there exists an identity 
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The K and 11 functions with this specific arrangement of the arguments, odd 
as it may look, are not unfamiliar to particle physicists. In fact these func- 
tions have been encountered in the analytical expressions for the phase space 
of three equal-mass particles [SI EJ [7] . To be precise, it is worth pointing out 
that the coefficient of 11 function accidentally vanishes in this example. Note 
the 3-body phase space can be obtained from cutting the simplest scalar two- 
loop sunset diagram. For a general sunset diagram with a nontrivial vertex 
structure, or if the power of propagators exceeds than 1, the complete ellip- 
tic integral of the third kind will inevitably arise in evaluating its imaginary 
part. Equation (2) can then be invoked to trade the 11 function for the K 
function, thus considerably simplifying the answer. 

In the aforementioned physical application, simple kinematics enforces 
< X < |, for which both the parameter and characteristic of K and 11 are 
positive and less than 1. This restriction is not necessary for general purpose, 
therefore it will be discarded in the following discussion. 

There exist some formulas which transform one 11 function to another 11 
function plus a K function [H |2]. However, the relation designated in (2), 
which links one 11 function to one K function only, is rather peculiar. To the 
best of my knowledge, this relation cannot be derived by any known formula, 
and it also has never been explicitly stated in any published work. For this 
reason, I feel it may be worthwhile to report it here. 

The strategy of the proof is to construct a differential equation satisfied 
by the left side of equation (2), called y{x) in shorthand. Employing the 
well-known differential properties of complete elliptic integrals [H H] : 
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together with the chain rule, I find that y satisfies the following first-order 
differential equation: 



, 1 + 22; + 3x2 
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The magic is that, after the differentiation, the complete elliptic integral of 
the second kind cancels, and the elliptic integrals of the first and the third 
kind conspire, in a rather peculiar way, to cluster into the original form. 

This is a basic type of linear differential equation, and the corresponding 
solution is 

y{x) = C , (5) 

X 

where C is a constant to be determined. Notice the above solution possesses 
a pole at X = and two branch points located at x = 1 and x = — |, and in 
general one should not expect C will assume a universal value in the entire 
domain of x, R. I shall attempt to fix the value of C region by region. 

First let us consider the case when x belongs to the open interval /q = 
(0, 1). Inspecting Eq.(2), it is easy to see y(xo) = for Xq = ^ G Iq. This 
initial value cannot be satisfied unless if C = 0. Also note the right hand side 
of (4) is a continuous function of x in this interval. Therefore, by the existence 
and uniqueness theorem for linear equation (for instance, see theorem 2.1 in 
Ref. |8]), the differential equation (4) admits the unique solution y = in 
this interval. 

Next I turn to the solution for x G Ji = {—00, — |)- Examining the left 
side of Eq.(2), one readily finds y{xi) = | for xi = — 1 G /i. This initial 
value can be satisfied only if C = — ^. Thus by the theorem of existence and 
uniqueness, the function in (5) with this value of C constitutes the unique 
solution in the interval /i. 

There still remain two other intervals, I2 = (— 1,0) and I3 = (l,oo) 
to be investigated. When x resides in I2, both K and 11 functions become 
complex-valued, and the left side of (2) turns to be purely imaginary; whereas 
as X G Is, though both K and 11 become also complex, the left side of (2) 
nevertheless is real. There are no simple initial values can be inferred in 
these regions. I performed a numerical check with the aid of the computing 
package Mathematica [4], and find the solution (5) with C = —j^ in these 
two regions agree with the left side of Eq. (2) to the 14 decimal place. 
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Equation (2) can have some interesting consequences. Here I illustrate 
one example: 
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The analytic expression for K{^) in term of F function is known [3J. However, 
it is amusing that these strange looking H functions can also be put in closed 
form. 

Lastly, one natural question may be raised- how about the validity of 
equation (2) when the domain of x is extended from real to complex? It is 
a curious question owning to the rich analytic structure of elliptic integrals. 
Through a numerical study using Mathematica, I find this relation still 
holds in most regions of complex plane. It is most lucid to demonstrate this 
examination in plots. The left side of equation (2) still vanishes in an oval 
region internally tangent to a rectangle (0 < Rez < 1, |Imz| < 0.33), which 
is characterized by the basin and plateau in Figure [H This can be viewed 
as a generalization to the first portion of equation (2). The second part of 
this relation is valid almost everywhere else except in a pear shaped region 
embedded inside a rectangle (1 < Rez < 6.74 and |Im2;| < 0.97), which 
can be clearly visualized in Figure [2] as the tip of an iceberg surrounded 
by the flat, boundless sea. At this stage I am unable to provide a rigorous 
justification for this observation. A thorough understanding of these features 
wiU be definitely desirable. 
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Figure 1: Profile of tlie left side of Equation (2) as a complex- valued function, 
with the real variable x promoted to a complex variable z. 
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